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ABSTRACT 

A new method of analyzing the 
kinematics of joint motion is 
developed in this study. Magnetic 
Resonance Imaging (MRI) offers 
several distinct advantages. Past 
methods of studying anatomic joint 
motion have usually centered on four 
approaches. These methods are x-ray 
projection, goniometric linkage 
analysis, sonic digitization, and 
landmark measurement of 
photogrammetry . Of these four, only 
x-ray is applicable for in vivo 
studies. The remaining three 
methods utilize other types of 
projections of inter-joint 
measurements, which can cause 
various types of error. MRI offers 
accuracy in measurement due to its 
tomographic nature (as opposed to 
projection) without the problems 
associated with x-ray dosage. 

Once the data acquisition of MR 
images was complete, the images for 
this study were processed using a 3- 
D volume rendering workstation. In 
this study, the metacarpalphalangeal 
(MCP ) joint of the left index finger 
was selected and reconstructed into 
a three-dimensional graphic display. 
From the reconstructed volumetric 
images, measurements of the angles 
of movement of the applicable bones 
were obtained and processed by 
analyzing the screw motion of the 
MCP joint. Landmark positions were 
chosen at distinctive locations of 
the joint at fixed image threshold 


intensity levels to ensure 
repeatability. The primarily two- 
dimensional planar motion of this 
joint was then studied using a 
method of constructing coordinate 
systems using three (or more) 
points. A transformation matrix 
based on a world coordinate system 
described the location and 
orientation of a local target 
coordinate system. 

This study will lead to future 
research involving volume rendering 
of MRI data focusing on the internal 
kinematics of the hand's individual 
ligaments, cartilage, tendons, etc. 
Its findings will show the 
applicability of MRI to joint 
kinematics for gaining further 
knowledge of the hand-glove (power- 
assisted) design for EVA. 

INTRODUCTION 

The purpose of this research is 
to determine the feasibility of 
using Magnetic Resonance Imaging 
(MRI) as a tool for the study of 
hand kinematics. This project was 
funded by the NASA Office of 
Aeronautics and Space Exploration 
and the NASA National Space Grant 
College and Fellowship Program. 
Various methods of determining the 
internal motion of the hand's 
internal bone structures exist, but 
MRI offers several advantages which 
will be discussed in detail later in 
this article. 
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DATA ACQUISITION METHODS 

The data for this study came 
from the images of the left hand of 
a 29 year old male. The entire hand 
was imaged , but for the purposes of 
this study, only the index finger 
was analyzed. The study focused on 
the rotation and translation 
characteristics of the proximal 
phalanx and its relation to the MCP 
joint, so that only the portions of 
the images pertaining to these areas 
were used. A restraint device was 
constructed to ensure that both 
random hand motion and vibratory 
noise were minimized during imaging. 
Although random hand motion can vary 
from one individual to the next, the 
magnitude of this motion is very 
small. Over the image acquisition 
period, this type of motion results 
in negligible position error. The 
device was designed to ensure that 
the index finger was separated from 
the other fingers in a fixed 
position. Six preset angles of 
rotation were possible using this 
device, of which three were used for 
this study. Since MR images' signal 
intensity is based on density of 
certain >atoms (hydrogen in this 
case) , markers containing water were 
placed on the restraint. Two of the 
markers were attached to the bottom 
of the device, while the third was 
attached to the hinged portion, just 
past the tip of the index finger 
(Fig 2) The markers appear as 

circles in Figure 2, with A denoting 
the marker beneath the hinge, B 
showing the marker affixed to the 
tip of the base, and C showing the 
marker at the fingertip. These 

markers allowed calibration of the 
procedure and a means for 
determining the accuracy of the 
group's results. 

MRI is used widely for its 
tomographic capabilities, high 
spatial resolution, and ability to 
produce high contrast between 

different tissue types. The physics 
of MRI are among the most 

complicated of all the current 
imaging modalities [2,3,7]. 
Briefly, the basic MR process is to 
place a subject into a strong 


uniform magnetic field in the bore 
of a large magnet, which typically 
has a field strength on the order of 
1.0- 1.5 Tesla. This field aligns 
the net magnetic spin of the 
individual atoms. In order to 
produce a usable signal, these 
aligned magnetic spins must be 
perturbed. A radiofrequency (RF) 
pulse is applied at the resonant 
frequency of the atomic spin of 
interest (hydrogen is usually chosen 
due to its abundance) . The spins 
are moved away from their preferred 
axis of rotation by the energy 
imparted by the RF pulse. The spins 
will try to realign with their 
original axis of rotation (along the 
bore of the magnet), and as they 
realign, they emit a RF signal. 
This signal can be decomposed into 
its frequency and phase components 
which are related to the location in 
the subject from which the signal 
was measured. The magnitude of the 
signal is proportional to the signal 
intensity generated by the various 
structures in the image. This 
magnitude is, in turn, related to 
the density of hydrogen atoms in the 
tissue and the realignment 
(relaxation) parameters of the 
tissue studied. These parameters 
are called Tl, the longitudinal 
relaxation time , and T2 , the 
transverse relaxation time. The 
density of the hydrogen atoms 
determines the material present in a 
give 3— D volume of space. This 
information is then projected onto a 
2-D display for viewing (Figure 1). 
The image data can be manipulated to 
reorient the slices in different 
orientations for further analysis. 

The programs used to control 
the MR scanner to create images are 
called pulse sequences. Each type 
of pulse sequence provides different 
contrast between various tissue 
types or has different spatial 
resolution or signal-to-noise 
properties. The most traditional 
(and slowest) type of pulse sequence 
is the spin echo pulse sequence. 
There are also many pulse sequences 
designed to provide essentially the 
same tissue contrast or signal-to- 
noise ratio as spin echo sequences, 
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3-0 Volun# 2-0 R«pr Mentation 



Figure 1 - Diagram showing the 3-D 
slices of a volume depicted on 
a 2-D CRT display. 


but in a fraction of the time. 
Since our goal was to optimize 
speed, image quality, resolution, 
and tissue contrast, we used a 
variation of a fast pulse sequence 
called TurboFLASH. The scanner used 
was a 1.5T clinical unit (Siemens 
Magnetom SP, I selin, NJ) . In order 
to obtain approximately isotropic 
image volume elements (voxels), we 
used a volumetric acquisition 
version of TurboFLASH known as 
MPRAGE (Magnetization Prepared Rapid 
Acquisition Gradient Echo) . With 
this sequence, we obtained 
acceptable spatial resolution 
( 1-2 mm) images through the entire 
hand with good differentiation 
between tissue types, and an 
acceptable signal-to-noise ratio in 
approximately 7 minutes . The pulse 
sequence parameters used were TR = 6 
msec, flip angle ■ 10 degrees, 

acquisition matrix = 130 x 256 x 
128, and a field of view = 250 mm. 

Figure 1 demonstrates the 

tomographic capability of MRI . 

Since a volume of image data of up 
to 30 cm 3 can be acquired using the 
particular radiofrequency coil 

chosen for this study, imaging the 
hand in an EVA glove will be a 
simple extension to this study. 


3-D VOLUME PREPARATION 

After acquiring the MRI data, 
it was necessary to prepare the 
volumes. These were created by 
adding the area of interest (the 
phalangeal bones of the index 
finger) for each slice of MRI data. 
Since each slice contains a 2-D 
representation of the phalanges, 
adding them creates the composite 
3-D volume. All volume rendering 
was completed utilizing a DMI 
MAXIVIEW 3-D dedicated workstation 
(MAXIVIEW, Dimensional Medicine, 
Inc., Minnetonka, MN) and its 
accompanying software. 

The MAXIVIEW workstation was 
created specifically as a tool for 
processing medical imaging modality 
data. It accepts image data from 
Computed Tomographic (CT) systems, 
Digital Subtraction Angiography 
(DSA) systems, and MRI systems. 
This system allows 3-D reformatting, 
resectioning and multi-planar 
reformatting for processing images. 
In addition to creating the volumes, 
MAXIVIEW also proved valuable for 
viewing, calibrating, and measuring 
the volumes. 

The overall problem with the 
thresholding of MRI data during 
volume rendering deserves further 
discussion. The reconstruction 
sequence of computer workstations 
such as the MAXIVIEW depends largely 
upon its ability to extract as much 
surface detail as possible in two 
dimensions (ie, differentiate 
between the surface of the proximal 
phalanx and the surrounding tissues 
in a given slice) and to then add 
these surfaces for a 3-D 
representation (Vannier and Geist, 
[10]) using MRI, CT, or DSA data. 
However, as the resolution of the 
individual slices decreases, the 
resolution of the corresponding 
volume likewise decreases. This 
fact makes CT an excellent choice 
for rendering. It is based on 
attenuation values for X-rays, 
called Hounsfield numbers, which 
offer a stark contrast between bone 
and soft tissue. However, CT still 
depends on X-ray attenuation, which 
makes MRI attractive for repetitive 
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studies of this type, owing to its 
lack of ionizing radiation. But MRI 
suffers from its lower resolution of 
the border between bone and soft 
tissue (compared to CT) , which often 
leads to difficulties in 
thresholding. This will become less 
of a problem as the contrast of MRI 
scanners' images increases. 

Thresholding the images is the 
key step in rendering volumes. The 
images contain gray scale computed 
Information based on a combination 
of hydrogen density, Tl and T2 
relaxation times (MRI) or Hounsfield 
numbers (CT) . This gray scale 
information is then utilized to 
determine what material is present 
inside a 3-D voxel. Voxels which 
contain the correct values are 
retained for inclusion in the 
volume, while others remain blank. 
It follows that the size of the 
voxel itself proves to be an 
overall limit to the volume's 
resolution. In this study, the 
typical volume was approximately 0.4 
nun in each dimension (x,y,z). 

After the volume is computed, the 
next step is the shading process. 
An initial viewing direction is 
chosen. ,* The 3-D view in this 
direction’ is projected onto a 2-D 
plane, called the image plane, for 
display on the CRT. 

Rendering the phalanges for 
this study was the next step. In an 
effort to maximize the resolution of 
the volumes, the markers and the 
bones were rendered separately. 
This created a separate 
representation for each, but the 
positions relative to each other 
remained fixed in the global 
coordinate system of the MRI 
scanner. For each slice of a pose, 
the appropriate bones were selected 
and thresholded. A second set of 
bone markers was added to the bone ' s 
rendered volume to ensure that the 
volume was correctly aligned while 
measuring in the 3-D mode (see 
Figures 3-5). After the volumes for 
the bones were created, the volumes 
for the restraint device's water 
markers were computed. The water 
markers served to calibrate and 
validate the measurements and 


algorithms for this study. A second 
set of markers were also produced 
with the water marker volumes to 
ensure alignment, just as we did 
with the phalanx volumes. 

In order to analyze the motion 
of rigid bodies, we used the 
calculations and methods devised by 
several authors [4,8,9] . The 

spatial coordinates of landmark 
positions (a minimum of three on 
each rigid body in each position) 
can be calculated from the 
measurement tools of the DMI. These 
coordinates then comprise the 
information necessary to derive the 
transformation matrix describing the 
motion, which is taken to be a screw 
motion. A screw motion can be 
described by a rotation matrix and a 
translation vector, or by the 
position of the screw axis, the 
angle of rotation about this axis , 
and the translation along this axis. 
This method is a fairly robust 
method of determining the motion of 
the rigid body, but does suffer from 
errors induced by scaling errors and 
random measurement errors (noise) . 

While a comprehensive 
treatment of the derivation of Spoor 
and Veldpaus' methods is not 
feasible in this study, several key 
points need to be clarified. 
First, for this method, three or 
more non— collinear points are needed 
in the first and second positions, 
with no initial approximation 
needed. The translation vector (v) 
and the rotation matrix (R) describe 
the movement of a rigid body from 
position 1 to position 2. The 
rotation matrix has the following 
properties : 


where I is the 3x3 identity matrix 
and T denotes transposition. The 
vectors describing each point in the 
first position (a i ) and those 
describing the points in the second 
position (q i ) are related by 

g i =Ra i +v (2) 

Since the unknown quantities in (2) 
are R and v, a method must be used 
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Fig 2 - A representative image of a "slice" 
of the left index finger in the restraint 
device. A, B, and C show the markers. 
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POSE #1 - COMPUTED UOLUME 



Fig 3 - The rendered volume in the first pose. 
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Fig 4 - The rendered volume in the second pose. 


POSE #3 — COMPUTED UOLUME ^ 



Fig 5 - The rendered volume in the third pose. 
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to obtain values for them by using 
the measurements of the points in 
the first and second positions. A 
computer program was written 
(Peterson, [8]) which takes the 
three-dimensional coordinates from 
two positions of the rigid body and 
then computes the transformation 
matrix in the form of 


|*ii *ii *ii °| 


T- 


*21 *22 *23 
*21 *32 *23 

Vj Vj v k 


0 

0 

1 


(3) 


jb/=nvax ( b*b l , b?b 2 , bfby ) ;n=± 


ft 7 *! 


(7) 


where n will be in the same 
direction as that derived from (4), 
but a scaling factor different in 
magnitude. The values for a point 
on the axis and the translation on 
the axis are then found from 


t =n T *v 


( 8 ) 


where the 3x3 matrix in the upper 
left describes the rotation (R) of 
the body, and v denotes the 
translation vector. Once the 
transformation matrix is determined, 
the methods of analyzing the motion 
can be one of several in existence. 
The method of Spoor and Veldpaus 
utilizes the symmetrical and anti- 
symmetrical portions of the rotation 
matrix to determine the body's 
motion about the helical axis, with 


and 

»=(-l/2).n.(n.v)*_ 2i2l|l 5 _. n .v (9) 

where s represents a radius vector 
of a point on the helical axis, and 
t is the translation along the 
helical axis. Noteworthy is the 
fact that the translation vector is 
quite sensitive to measurement 
errors and noise, while the rotation 
matrix is rather robust and yielded 
highly accurate results. 


2*n*sin(0) 




(4) 


The angular rotation is given by 

sin (0) -1/2 ,/ (Kjj-Rjj) J ♦ (K u --*ii> <* ji-*u> J ( 5 ) 


The rotation is also given by 

cos (6)- (1/2) ♦/?„-!) (5a) 


DATA ANALYSIS /RESULTS 

The angles of rotation about 
the helical axis for the respective 
position changes for this study were 


ANGLES OF ROTATION - POSE #3 
MARKERS 

Restraint Device - 52.0 deg 


Now that the angle is known, 
returning to equation (4) yields 
(n) , the unit vector along the 
helical axis. 

For the translation analysis, 
one must again derive an expression 
based upon the transformation 
matrix. The translation vector, v, 
is taken from the transformation 
matrix. The following information 
determines the necessary 
information: 

[l-cos (0) ] nn T - [ (R+R T ) -cos (0) I] /2 i g \ 

where bj, are three element columns. 
Another value for n is determined 


2- D Slice Image - 52.9 deg 

3- D Volume - 52.7 deg 


Table 1 - Calibration Results 


quite accurate, as were the 
calibration results utilizing the 
markers. Each marker set for the 
poses was rendered into a 3-D 
volume, just as the phalanx voxels 
were rendered. However, the water 
marker volumes were created for one 
slice of image data only. This was 
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to ensure that a clear rotation 
about one axis alone was visible, 
thus further validating the 
measurements of the DMI and the 
accuracy of the transformation 
matrix software. From the slice 
containing the markers ( see Figure 
2), measurements of the markers for 
pose #3 were determined for the 
angle B-A-C. This angle was then 
compared to the angle produced by 
measuring landmark positions from 
the rigid body A-B (notionally 
position 1) moving to A-C 
(notionally to position 2). The 
results are listed in Table 1. 
Furthermore, the unit vector along 
the helical axis for the rotation of 
Pose #3 was n = [0 0 1] T , thus 

showing that for the markers, a 
simple 2-D rotation occurred in the 
z-axis, which proves valuable later. 

After checking the process 
using the water markers, we then 
measured and analyzed the rigid body 
motion of the landmark positions on 
the proximal phalanx. The 

measurements were taken in the 
(x,y,z) coordinates and then 
transformed into the transformation 
matrix. From this, as in the case 
of the markers, we determined the 
motion of the phalanx about the 
helical axis using Spoor and 
Veldpaus algorithms. Seven landmark 
positions defined each rigid body 
motion , with future efforts geared 
toward increasing the number of 
landmark positions for increased 
accuracy and reliability. The 

results are listed in Table 2. 


From-To 

(deg) 

1-2 
2 - 3 
1-3 


12 3 


21.3 20.4 21.3 

12.4 14.0 12.7 

32.6 32.9 33.5 


Table 2 - The 3-D rotations 
for three trials. 


The overall consistency of the 
results above can be readily 


observed. The fact that the accuracy 
of the transformation increases as 
the magnitude of motion increases is 
a point of interest for future work. 
Also, as in the case of the markers, 
the motion was primarily two- 
dimensional about the z-axis, as 
expected. The small rotation in the 
x and y axes can be attributed 
mainly to the alignment of the 
subject hand. Random measurement 
errors also played a role . The 
mean and standard deviation of the 
movements are listed in Table 3. 


From-To Mean Deviation ( deg ) 

1-2 21.1 -8 

2-3 13.0 .7 

1-3 33.5 .4 


Table 3 - Statistical 
Analysis. 


CONCLUSION 

As noted earlier , future work 
in this area includes ongoing power- 
assisted glove work. A subject 

could wear the EVA glove inside the 
magnet, placing the hand within the 
radiofrequency coil. By stopping a 
subject at various stages of the 
motion required to perform a task, 
high resolution imaging of the 
entire had and glove volume could be 
performed at each stage. The 

transformations of the entire hand s 
members could then be analyzed. The 
effect of using the glove compared 
to non-glove motion could also be 
studied. Other areas of possible 
interest include those involving 
potential applications of injury 
related analysis of the hand. One 
major subject in this field is study 
of individual joint kinematics 
leading to injury suffered during 
certain repetitive activities, such 
as those encountered by a computer 
programmer while entering data at a 

keyboard . . 

Our research proved that MRI is 
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a feasible tool for kinematic 
analysis of rigid bodies. The 
process was very slow for this 
initial study, due to several 
factors, but subsequent studies will 
be completed much more quickly. The 
screw motion analysis of the 
proximal phalanx in this research 
produced results that were both 
reliable and repeatable for the 
proximal phalanx and should yield 
the same results for whole hand 
analysis, a future goal of the 
group. While other methods do 
exist, such as electrical 
goniometry, computed tomography, 
etc., only MRI offers a safe means 
of producing accurate results for in 
vivo analy s i s . 
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